Light scattering by ultracold atoms in an optical lattice 
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We investigate theoretically light scattering of photons by ultracold atoms in an optical lattice in 
the hnear regime. A full quantum theory for the atom-photon interactions is developed as a function 
of the atomic state in the lattice along the Mott-insulator - superfluid phase transition, and the 
photonic scattering cross section is evaluated as a function of the energy and of the direction of 
emission. The predictions of this theory are compared with the theoretical results of a recent work 
on Bragg scattering in time-of-flight measurements [A.M. Rey, et at, Phys. Rev. A 72, 023407 
(2005)]. We show that, when performing Bragg spectroscopy with light scattering, the photon recoil 
gives rise to an additional atomic site to site hopping, which can interfere with ordinary tunneling 
of matter waves and can significantly affect the photonic scattering cross section. 



I. INTRODUCTION 

Bragg scattering in condensed matter is a powerful 
method for gaining information over the structural prop- 
erties of crystalline solids. Usually, one employs ther- 
mal neutron beams, whose thermal wavelength is of the 
order of the interparticle distance inside the crystal. 
While elastic scattering allows one to measure the recip- 
rocal lattice primitive cell, inelastic scattering gives in- 
formation about the phonon spectrum and anharmonic- 
ities ^J. In atomic systems, Bragg scattering has been 
applied for demonstrating long-range order in structures 
of cold ions in traps and neutral atoms in optical lat- 
tices J,, 4,, _5, ^, 0] ■ Moreover it has proven to be a precise 
tool for the measurement of the elementary excitations 
of trapped Bose- Einstein condensate M, [§| and strongly- 
correlated atoms in optical lattices [13]. The spectra of 
the scattered photons, moreover, provide information on 
the details of atom-photon interactions. Studies on opto- 
mechanical systems, for instance, showed that the Stokes 
and anti-Stokes components of the scattered light may 
exhibit entanglement, which emerges from and is me- 
diated by the interaction with the quantum vibrational 
modes of the scattering system JJjJ . Such correlations are 
endorsed by quantum interference in the processes lead- 
ing to photon scattering, which is mainly visible in the 
height of the spectral peaks as a function of the emission 
angle , and can be an important resource for quantum 
networks 

In this paper we investigate the opto-mechanical prop- 
erties of strongly-correlated atoms in optical lattices. 
These systems present peculiar features, when compared 
with solid-state crystals. In optical lattices the bulk peri- 
odicity is determined by the light potentials, and is hence 
of the order of half the laser wavelength. One remark- 
able property is that light both couples to the atomic 
transition and is diffracted by the crystalline structure 
which the atoms form 0, |3j This property implies, 
for instance, that the system may exhibit peculiar self- 
organization, being the atoms a diffracting medium for 
the light which traps them 0, [H, [11, [l7| . 



In the dispersive regime, when the optical lattice can 
be considered a conservative potential, various states of 
ultracold matter can be realized [l^, thereby mimick- 
ing solid state models (lol . [20j , a prominent example of 
which is the quantum phase transition between a Mott- 
insulator and a superfluid state [2lj . Bragg spectroscopy 
provides an important tool for characterizing the quan- 
tum state of the atomic gas 8, 9, lOj. The experimental 
procedure typically uses two laser beams, whose wave 
vector difference q gives, by means of the mechanical ef- 
fects induced by photon recoil, a momentum and energy 
transfer hq and fku Q. The corresponding atomic re- 
sponse is detected by a time-of-flight measurement, con- 
sisting in releasing the trap and measuring the momen- 
tum distribution by atom detection An alterna- 
tive procedure makes use of parametric amplification fol- 
lowed by time-of-flight measurement, thereby revealing 
the energy transfer and the spectrum p^ . [2^ ]. These 
procedures mayallow one to measure the structure form 
factor [U, H^, and characterize the state of the gas. 

Most recently, ultracold atoms were loaded inside of 
optical resonators, and first measurements of the spec- 
trum of transmission of the light at the cavity out- 
put showed novel features, which can be brought back 
to the collective and coherent interaction of the atoms 
with the hght [13, [H [H, [H, [sO]. Several theoreti- 
cal works pointed out that the observation of the pho- 
ton scattered by ultracold atoms may provide com- 
plementar y in formation on the quantum state of the 
atoms (3Tl.l3^l33^33^35j. which could be non-destructive 
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in some setups 

We also remark that theoretical studies on Fermionic 
systems in an optical lattice showed that intensity fluc- 
tuations of the scattered light may allow one to deter- 
mine the temperature of the atomic cloud [37|. Optical 
detection, and in particular the intensity of the Bragg 
peaks, were proposed as a mean for revealing fractional 
particle numbers of Fermi gases confined by optical lat- 
tices [11]. In this paper, we study light scattering by ul- 
tracold bosonic atoms in an optical lattice, in the setup 
sketched in Fig. [T] We use a full quantum description of 
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FIG. 1: Light scattering by atoms trapped in a one- 
dimensional optical lattice with lattice constant do- The 
atoms are probed by a laser beam, with wave vector and 
frequency ujl, which couples to the atomic dipole transition 
at frequency cjq with ground state \g) and excited state |e) 
(see inset). The spectrum of the scattered light is measured 
at a detector as a function of the angle of emission. In exper- 
iments, one can also use a second laser beam, into which the 
photon is emitted with high probability, hence implementing 
stimulated Bragg scattering 



the photonic and atomic fields, for a range of optical lat- 
tice depths which covers the superfluid to Mott-insulator 
transition. By starting from the general Hamiltonian, we 
carry out the tight-binding and single-band approxima- 
tions, and we determine the scattering cross section of 
photons in the linear response regime. Extending previ- 
ous works [121, we systematically take into account the 
finite tunneling rate in evaluating the scattering cross 
section for parameters sweeping along the phase transi- 
tion Mott-insulator to superfluid state. Our study fo- 
cusses on a small lattice of 7 sites, and solves numeri- 
cally the Bose-Hubbard model for this system. In order 
to get insight into the numerical results, we also develop 
an analytical theory, which extends the theory presented 
in [23] by including the hopping induced by photon re- 
coil. The interference between the finite tunneling rate 
and the photon-induced hopping is visible in the height 
of the Stokes peaks as a function of the emission angle 
and can be revealed experimentally. 

This article is organized as follows: In Sec. [ll] we 
present the theoretical model. In Sec. IIIII the scattering 
cross section is evaluated both analytically and by means 
of numerical simulations. The conclusions are discussed 
in Sec. ITVl 



dom of the atoms are the electronic ground state \g) and 
an excited state |e) that form a dipolar transition with 
dipole moment D at the optical frequency ujq, which cou- 
ples to a weak laser probe. The Hamiltonian in second 
quantization reads H = 7?at + ^^omf + ^^int with [3^] 

i/at ^hiuoj dr^|(r)V'e(r) + ^ + H,g , (1) 

Hcnif = ^ hu;xa\ax , (2) 
A 

where by '0j(r) and i^ji^) we denoted the annihilation 
and creation operators of an atom in the internal state 
j = g,e at position r, and by a\ and the annihila- 
tion and creation operators of a photon in the mode at 
frequency lo\, wave vector and polarization eA -L k^. 
The atomic field operators obey the bosonic commuta- 
tion relations [iljj{r),7pjr{r')] = [■(/;j(r), ■(/'j, (r')] = and 
[V'j (r), -0^, (r')] = Sjjr5{r—r'). The Hamiltonian term 
(He) describes the motion of the atoms in the internal 
state \g) (|e)), and H^g gives the coUisional interaction 
between the atoms in states \g) and |e). We will assume 
that the atoms interact with radiation far-off resonance 
from the dipolar transition, hence the occupation of the 
excited state is small and will be neglected. Therefore, 
we just need to provide the detailed form of the ground 
state term. 
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^ I dryjlir)^lir)ijgir)^gir), (3) 



where Ugg is the strength of the contact interaction. The 
potential V{r) is assumed to be periodic along the x- 
direction and reads 



V{r) = Vo sin^ j + ^muj^iy^ + z^) , (4) 

where Vq is the lattice depth, do the lattice constant, 
and ujr the frequency of the harmonic trap which tightly 
confines the transverse motion. 

Finally, the interaction term between atoms and light 
reads (in the length gauge) 

i?int = Y.^^>^j dr^I (r)V<,(r)aAe*-'- + H.c. , (5) 



where 



II. THE MODEL 

The scattering system we consider is composed by N 
identical bosonic atoms of mass to in a periodic potential, 
as shown in Fig. [1] The relevant internal degrees of free- 
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is the coupling strength, with eg the vacuum electric per- 
mittivity and V the quantization volume. 
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A. Linear response 

At room temperature and equilibrium, the atoms are 
in the electronic ground state and the state of the optical 
modes of the electromagnetic field can be approximated 
with the vacuum |0). We now assume that a laser, at 
frequency wl and wave vector Icl , couples to the atomic 
dipole transition. The laser field is described by a co- 
herent state of the corresponding electromagnetic field 
mode with amplitude a^, such that the mean number 
of photons is given by |aLp- In the regime in which the 
atom-laser coupling is sufficiently weak, corresponding to 
the condition ICiaij <C |cjo — we eliminate the ex- 
cited state from the equations of motion of the ground 
state in second-order perturbation theory in the small pa- 
rameter |Cj;,ai|/|a;o — w^l- The dynamics of the atoms 
in the electronic ground state \g) is now described by the 
effective Hamiltonian 



where we denoted by |1a) = a>,|0) the state of the elec- 
tromagnetic field with one photon in mode A, and by \i) 
and I/) the states of the atoms before and after the scat- 
tering, respectively, which are eigenstates of Hamiltonian 
Hg at energies Ei and Ef. The function 

sin(a;r/2) 

is the diffraction function, giving energy conservation for 
infinite interaction times, limr^oo 5'^^\uj) — 5{lo) [43j. 

Equation (fTO|) shows clearly that the scattering rate de- 
pends on the state of the atoms before and after the scat- 
tering event. In the following, we derive the atom-light 
interaction Hamiltonian in the tight-binding approxima- 
tion and conclude this section by introducing the many- 
body atomic states which are relevant for the scattering 
process considered here. 



where the interaction term takes the form 
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ax' / dre"i"-^/.J(r)Vg(r) 



etc. 



A, A 



(8) 



and describes the absorption of a photon in the mode 
A' and wave vector k^'and the emission into the mode A 
and wave vector k^, weighted by the Fourier transform 
of the density AAq = / dre"i '"?/'|(r)'!/'g(r), with 



q = kA' - k^ 



(9) 



In the following we will assume that the interaction be- 
tween photons and atoms is essentially Hamiltonian, and 
hence fully determined by the Schrodinger equation gov- 
erned by Eq. ([7])- This is valid in the regime which we 
consider in this article, namely, when the detuning of the 
light \ll}q — ujl\ S> 7, with 7 the linewidth of the excited 
state. 

In this work, we study Bragg scattering of laser pho- 
tons by atoms in the one-dimensional periodic array given 
by potential ([4]). We will hence evaluate the differential 
scattering cross section for coherent scattering. Assum- 
ing that jail ^ 1, so that the atoms absorb at most one 
photon from the laser at a time, the differential scatter- 
ing cross section is found from the rate of scattering one 
laser photon into the mode A. In particular, the scatter- 
ing rate reads 



27r 



Tx^^X = |(lA,/|i?L|lAi,,«)| 
/ 

X d^^\LUL~u;x + iE,~Ef)/h), 



(10) 



B. Tight-binding regime 

We assume that the atomic wavefunctions are well lo- 
calized at the lattice minima, such that the tight-binding 
approximation can be applied. Furthermore, at ultralow 
temperature and not too strong interactions, the atomic 
gas is in the lowest band of the periodic potential and 
in the ground state of the radial oscillator, so that the 
atomic field operator can be decomposed as 



(12) 



where wi{x) = w{x — 
tered at position Wq, 



Mq) is the Wannier function cen- 
with the sum going over all lat- 



tice sites, and (j)o{p) = exp{—p'^/2^^)/{^rV^) is the 
ground state of the radial oscillator (p — y^y^ + z^) with 

— y^h/mojr. The operators bi annihilate an atom at 
site I and fulfill the standard bosonic commutation rela- 
tions = Si^i>. Using this decomposition in Eq. ([3]), 
allowing only nearest-neighbour hopping and restricting 
to on-site atom-atom interactions, we obtain the Bose- 
Hubbard Hamiltonian [131 



ni 



(13) 



where ni = b^bi is the atomic number operator at site I 
and /i is the chemical potential. The coefficients for the 
hopping term and the on-site interaction strength read 



J 



dxwiix) [ — + V{x) ] wi+i{x) , (14) 



U = WgalTF / <^xwi{x) 



(15) 
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with the Wannier functions chosen to be real. Note that 
frozen transverse dynamics, as assumed in Eq. Ijl2p . is 
here ensured by taking J, U{n) <^ hur, where (n) is the 
mean site occupation. Within this decomposition the 
term describing atom-hght scattering takes the form 



H', 



E 



A, A 



(16) 



Here, 



r(q) ^ ^e'«^'''« [jo(q)n, + Ji(q) + bj^.bi 

(17) 

consists of a photon-dependent energy shift, weighted by 
the coefficient 



Jo(q) = e-3('?'+9')«' y da;e''?="^^^;o(a;)^ (18) 
and a hopping term with coefficient 

Ji(q) = 6-3(9'+'?')?' y da:wo(a:)e'«-^'u;o(a;-do)(19) 

which describes hght-assisted tunnehng due to the me- 
chanical effects of photon scattering. This latter term has 
been neglected in previous theoretical treatments [23 . [3^ . 
[33I I . Its effect has been investigated in Ref. [sil, lio] for 
light scattering by ultracold atoms in a double well po- 
tential, showing that the mechanical effect of light can in- 
terfere with ordinary tunneling between the wells, gener- 
ating observable effects in the first-order coherence prop- 
erties of the scattered light. We hence expect that it will 
give rise to observable effects in the Bragg signal by ul- 
tracold atoms in optical lattices. 

In the following, we introduce the many-body states, 
eigenstates of H'g, which are relevant for the scattering 
process when the system is in the Mott-insulator and 
in the superfluid regime. In this treatment we use the 
same notations as in Ref. [4l|, and refer the reader to 
this work for more details, like, e.g., the careful compari- 
son between the Bogoliubov approximation and the exact 
solution for small onc-dimcnsional systems. 

1. Mott-insulator State 

For vanishing hopping, the ground state of Hamilto- 
nian ([TH|) is the Mott-insulator state with all lattice sites 
equally occupied with (integer) filling factor g = N/M, 



(0)\ 



M 

n 

1=1 



|0)at-|5,5,-<?,5), (20) 



where |0)at denotes the vacuum. The corresponding 
ground state energy for J = is easily found and reads 



= MUg{g-l)/2-Mgn. The lowest-lying excitations 
take the form 



(0) 



V9{9 + 1) 



(21) 



where one particle and one hole are created at site n and 
TO, respectively, with energy — + U. These states 
form a degenerate subspace of dimension M{M—1). This 
degeneracy is lifted for finite values of the hopping J. 

The corrections due to a non- vanishing but small value 
of tunneling are evaluated using perturbation theory. In- 
cluding the first-order correction, the ground state now 
reads 



J2 

u 



jj^Mg(g + 1) 



(0)\ 

/ 



--^2Mg{g + l)\S). 



where \S) 



/2M 



(0) X 



(0) 



(22) 

_^)^ is the nor- 



malized state of adjacent particle-hole excitations, while 
the term at second order in J warrants normalization of 
state (P^ . The corresponding energy is Eq = Eq+0{,P). 
The lowest-lying excitations are determined using degen- 
erate perturbation theory within the subspace of single 
particle- hole excitations. 



(0) \ 
n.mi 1 



(23) 



where the coefficients c„Jm fulfill the normalization con- 
dition and satisfy the equations 



(5 + l)(ci+i_m + cJjL]^_„J -t- (7(cjj'„j^]^ -I- cjj^„j_]^) — Ai 

with periodic boundary conditions 



0. 



(24) 

(25) 
(26) 



Term Ai in Eq. ([M)) is the first-order correction to the 
corresponding energy, Ei — Eq + U ~ JAi -I- O(J^). 

An analytic solution of Eqs . (|24l) can be derived in the 
limit of large filling 5^1 [2J, |4l| . This limit introduces 
a symmetry between particle and hole excitations, that 
simplifies the analytical treatment but imposes a selec- 
tion rule, which is strictly correct only when g ^ 00. The 
coefficients, evaluated in this limit, read 



\/2 r sin [ar\n - m\] e'""("+™) 
M \ sin [ar{n ~ to)] e'""("+™) 



for r + s odd, 
for r + s even, 
(27) 

with a = s = 0, 1...M- 1 and r = l,2...Af- 1. 
Correspondingly, the lowest-lying excitations and their 
energy read (at first order in J and for g ^ 1) 



5 



(0), 



for s — and r odd, 



otherwise, 



E° + U- 2J{2g + 1) cos ar cos as + 0( J ) , 

I 



(28a) 
(28b) 



where Afr is a normahzation factor. Note that states 
contain a correction proportional to the ground 



state I ■00 



(0)\ 



This correction is found from non- 



degenerate perturbation theory and warrants the or- 
thonormahty of the new basis {\ipQ^''), \4'\r\])}- 



2. Superfluid state 

In the weakly-interacting superfluid regime, quantum 
fluctuations in the number of atoms per site are described 
by the decomposition 



Term H2 describes the dynamics of the non-condensed 
fraction at leading order. It can be written in the diago- 
nal form 



H2 = Y^ hrip 



1 



(33) 



where operators ap and are respectively the bosonic 
annihilation and creation operators of the Bogoliubov 
excitation with quasimomentum p — n2'K / Mdo, with 
n = —M, M ~ 1, . . . , A/ — 1, and the frequency Vlp is 
given by 



(29) 



where zi is a complex number describing the order pa- 
rameter and (3i is the fluctuations operator obeying the 
bosonic commutation rules. The order parameter zi is 
found by minimizing Hamiltonian IjlSp at zero-th order 
in the expansion in /3; and . It obeys the discrete non- 
linear Schrodinger equation 



JkZmSkJ + U\zi\'^zi , 



(30) 



(fc,' 



where \zi\'^ corresponds to the condensate fraction. For a 
translationally invariant lattice and in the limit of weak 
interactions, as considered here, it is given by = 
and Eq. ([50]) reduces to 



= -2J + Ug. 



(31) 



Using Eq. (|^^ and Eq. ([50)) in Hamiltonian , keeping 
only terms up to second order in the operators /3i , (3j , one 



finds H'g = Ho + H2 with 



H2=Y, ^Lm^lPm + Ml,rnpj PI + H.C. , (32) 

l,m 

where the coefficients read 

(n.k) 



with 



4Jsin^ 



(34) 



(35) 



Operators ap, alp satisfy the commutation relations 



Sp,p' and are related to j3i by the Bogoliubov 
transformation 



(3i = ^ Ve'f'''''w„a 



e-'P^o^-pat. (36) 



The Bogoliubov amplitudes Up, Vp satisfy the equation 
|upp — |fpp = 1, as a consequence of the commutation 
relations, and depend only on the modulus of the quasi- 



momentum. Up = U- 



pj 



-p. They are solutions of 



the Bogoliubov-de Gennes equations, which in our case 
read 



Jtp + Ug ~Ug 
Ug -j£p~Ug 



(37) 



In particular. 



Ep + Ug + hflp 

2hnp 

tp + Ug — hflp 



2hnr, 



UjiVr} 



2hn„ 



(38a) 
(38b) 
(38c) 
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We note that Cp is the energy of a non interacting par- 
ticle in the lattice. By replacing it with the free-space 
energy ep — s- /2m we recover in Eq. p4p the dispersion 
relation for a weakly-interacting dilute Bose gas in free 
space (42]. Contrarily to the case of the uniform one- 
dimensional system, where Bogoliubov theory is not ap- 
plicable, for a finite system it provides a well defined and 
small depletion for [// J <C 1 and large filling. The corre- 
sponding spectrum of the differential scattering cross sec- 
tion will be compared below with the numerical results 
obtained for a finite Bose-Hubbard model composed of 7 
atoms. 

In the following, we will denote by |0)sf the super- 
fluid state, where all atoms are in the condensate, and 
by |p)sF = ctp|0)sF the state with one Bogoliubov excita- 
tion at quasimomentum p. In particular, we will consider 
scattering processes, such that the state of the atoms will 
include at most one Bogoliubov excitation. To this aim, 
it is convenient to rewrite operator T(q) in Eq. (|17p using 
the decomposition of operator hi in Eq. 

rsF(q) = T^^ (q) -I- (q) + Tg> (q) , (39) 

The first term on the right-hand side of the equation de- 
scribes radiation coupling with the condensate and reads 

(q) = 9 ( Jo(q) + 2 Ji(q)) ^ e"^^"^o ^ (40) 

I 

while the other terms give radiation coupling with the 
Bogoliubov excitations, and take the form 

(q) = E e''^''" (^0 (/^' + ) (41) 
I 

+^i(q)(/3/ + A+i + /3/+i + /30) , 
^sF^q) = Y.'^'''-"'°{m^)P\Pi (42) 

I 

-HJi(q)(/3//3i+i+/3,ViA)) , 



where k = nfc and |/) are the initial and final atomic 
states, eigenstates of Hamiltonian (fT3|) at the eigen- 
frequencies uji and w/, respectively. Using Eq. ^ in 
Eq. one can easily verify that the differential scatter- 
ing cross section is proportional to the dynamic structure 
factor [H. 

We evaluate the scattering cross section assuming that 



where the superscript gives the order in the Bogoliubov 
expansion. 



III. LIGHT SCATTERING 



Light scattering by a one-dimensional optical lattice 
of ultracold atoms is studied in the setup sketched in 
Fig. [TJ A laser plane wave at wave vector kL , frequency 
lul = c|kL|, in a coherent state with amplitude ai,, drives 
the atoms. We evaluate the scattered light as a function 
of the angle of emission, determined by the wave vector 
k of the mode into which the photon is emitted, and of 
the frequency of the emitted photon. 

The scattering process is evaluated assuming that the 
laser very weakly excites the atom, so that the atom- 
photon interaction is described at lowest order by Hamil- 
tonian ([5]). More in detail, the condition ja^l <C 1 means 
that the atomic sample is driven by at most one pho- 
ton. A scattering process will then occur with proba- 
bility |q!lP and will consist of the absorption of one in- 
cident photon in the mode of the laser, represented by 
the state |1l), and the emission of a photon in one of 
the modes of the electromagnetic field at wave vector k 
and polarization Ck -L k, represented by the state |lk,e)- 
The corresponding differential scattering cross section for 
the photon scattered at frequency oj in direction n in the 
solid angle Vl is pro por tional to the scattering rate (flOl) 
and takes the form [43l . |4^ 



(43) 

I 

the atoms are initially in the ground state either of the 
Mott-insulator or of the superfluid phase, and that the 
atoms are scattered into a final state belonging to the 
lowest-lying atomic excitations. Using the form of oper- 
ator H[^^ in Eq. ([16]), Eq. (|43|) can be written as 



J 



'"("'^^ = jSA^ E E K/' ^^AHL\h <5(^)(c.i + _ ^ _ ^^) , 
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where 



.(1) 



(44) 



where G are the vectors of the (one-dimensional) recip- 
rocal lattice and 



.(Af) ^ 1 sin^(Af dog/2) 
" AP sm\doq/2) 



(49) 



a^''^ ^A{n)\{t\T{q)\t)\U^^\LOL-Lo) (45) 
gives the elastic component of the scattered light, while 

aW = ^(f])^|(/,zk|r(q)|j,li)|2 5(^)(^i-a;-,5a;/) 
/ 

(46) 

describes the scattering events in which one mechanical 
excitation at frequency Sujf is absorbed from the pho- 
ton by the atomic lattice (Stokes component) and corre- 
sponds to the one-phonon terms in neutron scattering [l| . 
The corresponding phonon emission processes, giving the 
anti-Stokes component, are here absent as initially the 
atoms are in the ground state. Moreover, higher order 
terms, corresponding to higher-order phonon terms in 
neutron scattering, are here neglected as we assume that 
at most one mechanical excitation is exchanged between 
lattice and photons. 

The operator T(q) in the above equations is given in 
Eq. P7|) . while the coefficient A{^1) depends on the angle 
of emission and takes the form 



gives conservation of the Bloch momentum in a finite 
lattice with M sites, such that S^^^ Sq^c (Kronecker 
delta) as M ^ oo. In Eq. (gS]) we omitted terms at third 
and higher order in J and Ji(q). This approximation 
will be applied to the rest of this section, assuming that 
these higher-order terms can be neglected. 

The presence of in Eq. expresses the von- 

Laue condition for Bragg scattering. At zero order in the 
hopping term, Eq. (I48p gives the response of a crystal 
of particles oscillating around their equilibrium position. 
In fact, using a Gaussian ansatz for the wave functions, 
one can estimate |Jo(q)P ~ e"^'^, with W = [g^^^ -I- 
(9y + 9z)fr]/8' where ^j. and are the widths the atomic 
wave functions in the axial and radial direction, showing 
explicitly that this term is analogous to the Debye- Waller 
factor P, Si]. The term proportional to J is instead 
a novel feature with respect to traditional condensed- 
matter systems, that arises from light induced tunneling. 

The Stokes component for the Mott-insulator is evalu- 
ated taking the final states |/) 
and reads 



iVf.'ij) given in Eq. ([28 



E 



1^ 
c A2 



Stt 



|D-n|^ 
IDP 



(47) 



aW(f7,c.) = Ain)J2sm'[^)\Br,s\' 



(50) 



X S{UJL -tJ-Wr,s)4ff),G' 



with LL}r,s = {Er,s ^ Eq) / H, and where we have introduced 



where 7 is the linewidth of the dipole transition, A ~ 
ujl — ^0 is the detuning of the laser from the atomic 
transition and Q,q = ^Juj^ l2heaTi ■ e^. 



Scattering cross section as a function of the 
atomic state 



We now give an analytic expression for the scattering 
cross section in Eq. ([43]) for the initial and final states 
determined in Sec. Ill B II and III B 21 



1. Mott-insulator 



For the Mott-insulator phase the initial state is \i) 
IV'o^^) given in Eq. Using Eq. (HH), we find 



a^^,{n,u;) A{m'S{iOL - c.)Cg (48) 



q{s) = qx 



27r 
Mdo * ■ 

The coefficient in Eq. (|50| takes the form 

Ji (q) for r H- s odd, 
2^Jo(q)sin(i^s 



(51) 



I Jo(q)P + 4V<?(5 + l)-Re{Jo*(q) Ji(q)} ) , 



for r + s even, 
(52) 

showing that the transition to the excited states with 
r + s odd is due to photon recoil, and is hence a light- 
induced hopping process. Note that condition q{s) = G 
shows that the quantum number s, and more specifically 
2tts/L, with L = Mdo the length of the lattice, plays 

the role of the quasi- momentum of the states \ipi]}). We 
remark that Eq. (|52|1. for r + s even, agrees with the 
result evaluated in [2J] (see Eq. (9) of that paper for 
comparison). The result we find for r -I- s odd, on the 
contrary, is discarded in the treatment of [13], as there 
the authors neglected light induced hopping terms. In 
the Mott-insulator regime these terms are usually very 
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small with respect to the other contributions. They give 
rise to a significant contribution when interfering with 
ordinary tunneling. This latter type of contributions is 
ruled out in the analytical model by the selection rule 
introduced by the assumption g ^ 1, but it is visible in 
the numerical results, as it will be shown in Sec. IIIII 



initially prepared in the Bose-Einstein condensate. In ad- 
dition, for the analytical calculation we consider the limit 
U 0. In this limit we can neglect the quantum deple- 
tion of the condensate due to the interactions and take 
the initial state \i) — |0)sf according to our notation. 
The zero-phonon term takes now the form 



2. Superfluid 

When evaluating the differential scattering cross sec- 
tion in the superfluid phase, we assume all atoms to be 



J 



^SF 



|Jo(q)+2Ji(q) 



2E- 



N 



-Re {( Jo(q) + 2Ji(q))*(Jo(q) + 2Ji(q) cos (pdo))} 



showing that the light-induced tunneling effects enter already at first order in this expression. As in the Mott-insulator 
case, the analogous of the Debye-Waller factor can be here identified in the term |Jo(q)p. In this case, though, 
tunneling effects become more important, modifying significantly the signal as we will show. The first-phonon term 
reads 



'SF 



f^p)^|(^o(q) + Ji(q)(l+e-'P*)) 



1 2 AM) 

I "q.-p,G ' 



r 



(54) 



and describes the creation of Bogoliubov excitations with 
quasi-momentum hp by photon scattering, such that the 
relation p = — G holds. 



B. Numerical results 

In this section we report the numerical results for the 
differential scattering cross section obtained when the 
atoms are in the Mott-insulator or in the superfluid state. 
The numerical results are obtained for a lattice of M = 7 
sites and fixed particle number N = M . The coefficient 
entering the Bose-Hubbard Hamiltonian in Eq. (fT3| and 
the operator T(q) in Eq. pT|) are calculated by using 
the Wannier functions relative to a given lattice depth 
Vb of optical potential Hamiltonian is diag- 

onalized exactly and the corresponding states are used 
for determining the differential scattering cross section in 
Eq. (|43p. The numerical results are also compared with 
the analytical predictions of the scattering cross sections 
reported in the previous section. Although the latter are 
valid for very large lattices and for large mean site occu- 
pation 5 3> 1, we find reasonable agreement when com- 
paring these predictions with those for a small lattice of 
7 sites and single occupancy (see also [2^1). 

Figure [2Ia) and (c) display the one-phonon contribu- 



tion to the differential scattering cross section, cr*^^^ (il, cj), 
as a function of the frequency lo and for different scat- 
tering angles when the atoms are in the Mott-insulator 
state. The numerical results are compared with the an- 
alytical model (dashed line) and with the model used in 
the numerical simulations in 24], in which light-induced 
hopping terms are not considered. 

We first discuss the numerical results which most 
closely approach the exact solution. The appearance 
of multiple peaks corresponds to the excitations of the 
atoms in the Mott-insulator due to the photon recoil. 
The number of peaks for the numerical result is M — I, 
which correspond in this case to 6. They can be individu- 
ally resolved, as the system considered here is finite, and 
the width of each individual peak is limited b y th e de- 
tection time T (or the spectral resolution 1 /T) [4^ . The 
analytical results are found using the model described in 
Sec. |IT1 which assumes a large on-site occupation. They 
are characterized by the same peak number, although 
only half of them is visible in the figure. In fact, the 
intensity of the peaks arising from the coupling of the 
ground state to the corresponding excitation via light- 
induced hopping (corresponding to the terms in Eq. (j52[) 
with r + s odd) are very small compared to the other 
ones (corresponding to the terms with r + s even) and are 
therefore not visible (note that, due to the assumption of 




q^do = 27r/7 qxdo = 




FIG. 2: (color online) Stokes component of the diflFerential scattering cross section (in units of A{Q,)) as a function of frequency 
(in units of the recoil frequency uir) for two different scattering angles, corresponding to Qxdo — 2-k jl (top row) and to 
</a;do = 67r/7 (bottom row). The curves have been evaluated for a lattice of M = 7 site and N = M = 7 composed by ^'^Rb 
atoms in the |F = 2, mp = 2) hyperfine ground state. The black solid line corresponds to the numerical results, the blue dashed 
line to the analytical formulas (see text), the red dashed-dotted line to the model of [2^ . where the light-induced hopping is 
neglected. Plots (a) and (c) are evaluated for Vb = 9>.lhujR {U / J ~ 17) which corresponds to the Mott-insulator state. Plots (b) 
and (d) are evaluated for Vb — O.lhuii ( (7/ J ~ 1) which corresponds to the superfluid state. Other parameters are do ~ 413nm, 
as = 105ao with ao being the Bohr radius, and lj,. = Wujr corresponding to the experimental parameters in .22] (For these 
parameters the size of the radial wavepacket is ~ lOos). The frequency resolution is set to Aw = 300 Hz, corresponding to 
an integration time T = 3 msec. 



large on-site occupation, interference between ordinary 
tunneling and light-induced hopping is suppressed). The 
central positions of the visible peaks present a systematic 
shift with respect to the ones found numerically. This 
systematic shift originates from the assumption g ^ 1, 



and has been observed in [24|. Nevertheless, the analyti- 
cal solution still provides some insight into the numerical 
results. In Eq. (pS)) . using Eq. (|?T|) we find that the peaks 
are centered around the energy E' = U with a spreading 
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FIG. 3; (color online) Stokes component of the differential scattering cross section (in units of A{^1)) as a function of the 
frequency (in units of ujr) and of the Bragg angle = qxdo (in units of tt). The plots have been evaluated numerically for (a) 
Vo = S.lhujR and U/J^ 17, (b) Vo — O.lfitJjj and U/J ~ 1, (c) Vo = O.lhuiR and U/J ~ 0.1. The other parameters are as in 
Fig. El 



about this mean value of width ^J{2g + l)cos (^^^^ 

Such spreading decreases as qxda approaches tt, compare 
Fig. HJa) and (c). In particular, for qxdf) — tt, the width 
of the distribution of the Stokes excitations vanishes and 
the spectrum reduces to a single peak, corresponding to 
the on-site energy U . 

The results for the superfluid regime are reported in 
Figs, mb) and (d). Here, the analytical solution predicts 
that in the limit g ^ \ the total momentum of pho- 
ton and lattice is conserved in a scattering event. Such 
property implies that the Bogoliubov mode matching the 
momentum-conservation condition, is excited, and there- 
fore one expects a single peak in the spectrum. For N = 7 
atoms and g — 1, the numerical results for U/J~l give 
a single peak at qxdo = 27r/7, while at qxdo = Gtt/T mul- 
tiple peaks are found. In this case, instead of a collective 
density fluctuation with a well defined momentum p, one 
observes particle- hole types of excitations as in the Mott- 



insulator case. In Fig. [2td) one observes a larger spread 
of the peaks as compared to the Mott-insulator case at 
the same Bragg angle. This is due to the larger value of 
the tunneling rate J. We remark that, choosing smaller 
values of the ratio U/J by ramping down the on-site in- 
teraction strength, as it is shown below, the spectrum re- 
duces to a single peak at all Bragg angles and approaches 
the limit of the single-particle spectrum, as it is recovered 
in Eq. §^ by setting f7 = 0. 

We now compare the numerical results, obtained tak- 
ing systematically into account the light-induced hop- 
ping term, to the results found when this term is ne- 
glected, corresponding to the treatment in 2j|. In the 
Mott-insulator case, comparison between the numerical 
results with and without light-induced hopping effects 
shows that in the first case one finds interference between 
ordinary tunneling and light-induced hopping. This gives 
rise to an alternating enhancement and reduction of the 
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q^do = 67r/7 




3 [0.3] 
Vo/Er [U/J] 



6 [0,9] 



9 [2,2] 



12 [4.7] 



15 [9,5] 



FIG. 4: (color online) Contour plot of the Stokes component 
of the differential scattering cross section (in units of A{0.)) as 
a function of the frequency (in units of ljr) and of the lattice 
depth Vb in units of u>r for qxdo = 6tt/7 (the corresponding 
value of the ratio U/ J is reported in the axis between squared 
bracket). The black dashed line marks the critical value at 
which the phase transition occurs in the thermodynamic limit. 
The other parameters are given in Fig. [2] 



peak heights at different frequencies, which is absent in 
the model discarding hght-induced hopping effects. In 
general, the light-hopping term contributes in determin- 
ing the height of some peaks, giving substantial modifica- 
tions of the spectrum which can be revealed experimen- 
tally. The effect is larger in the superfluid regime, where 
tunneling is enhanced, as one can see in Fig. [H^b). Here, 
the central peak at qxdo = 2t:/7 is 50% higher than in 
absence of this contribution. 

Figures [3ja)-(c) display the spectra of cr'^^^ as a func- 
tion of the frequency and of the Bragg angle, in three 
different points of the phase diagram. We remark that 
the width and spacing of the Bragg peaks are deter- 
mined by the finite size of the lattice. The plots in (a) 
and (b) are made in the same parameter regimes as in 
Fig. m (a),(c) and (b),(d), respectively, , namely U/J ~ 
17 and U/J ~ 1. Figure OJc), instead, corresponds to 
the value U/J^ 0.1. Here, one observes almost a single 
peak at each Bragg angle, as expected in the weakly- 
interacting superfluid phase. 

Figure HI shows c'^^-' as a function of the frequency 
and the depth of the potential, hence sweeping from 
the Mott-insulator to the superfluid regime at a given 
Bragg angle, corresponding to large momentum transfer 
{qxdo = 67r/7). Here, one observes that the spectrum 
varies from multiple peaks, deep in the Mott-insulator 
regime, to a single peak in the weakly-interacting super- 
fluid regime. The single peak appears around a value 
of U/J much smaller than the critical value [U/J]c for 
the Mott to superfluid transition (which, in the ther- 
modynamic limit, is predicted for [U/J]c = 3.37, see 
Ref. 47]). The presence of multiple peaks also in the 
superfluid phase close to the phase transition is remi- 
niscent of a strongly-interacting superfluid phase. Such 



phase contains, beyond the gapless phononic modes, also 
gapped modes [H, |4^ [s^, [5l|, [Hll , which are predicted 
to be dominant at large quasi-momentum. We expect 
that also in the thermodynamic limit the transition to 
a single peak in the scattered-light spectrum will occur 
at lower values oi U/J than the Mott-insulator to su- 
perfiuid phase transition and will be also dependent on 
the momentum transfer. The identification of the Mott- 
insulator to superfluid phase transition should rather rely 
on the existence of a gapless spectrum. In spite of the 
very small size of the considered system, indications of a 
gapless spectrum are present in our results, as one can 
see comparing Fig. [3l[a) with (b),(c). 

The intensity of the scattered light as a function of the 
Bragg angle is determined by the differential scattering 
cross section 



(55) 



and is reported in Figs. [5] for the atoms in (a) the Mott- 
insulator and in (b) the superfluid state. The solid line 
here corresponds to the numerical results, the dashed line 
to the analytical predictions and the dashed-dotted line 
to the model where light-induced hopping has been dis- 
carded, similar to the case considered in Ref. [l^]. In this 
latter work, in fact, corrections due to the tunneling J 
were neglected when evaluating light scattering by the 
atoms in the Mott-insulator state, while the calculation 
of light scattering from the superfluid state was made dis- 
carding the flnite value of the on-site interaction as well 
as the flnite width of the Wannicr functions. In Fig. [5l^a) 
we observe that in the Mott-insulator regime the signal is 
dominated by the elastic component, and corresponds to 
a classical diffraction grating. In the superfluid regime, 
on the other hand, one flnds that the amplitude of the 
Bragg peak is modified, and a background signal appears 
which is due to light scattering by the condensate frac- 
tion. This signal is the signature of the superfiuid phase, 
and it arises from the coherent effects of tunneling and 
light-induced hopping. We also notice that in the super- 
fluid phase only the flrst diffraction order is visible. This 
is due to the increased width of the atomic wavefunction, 
which yields a faster decaying Debye- Waller factor Jq (q) . 
We remark that higher diffraction orders would be visi- 
ble if the superfluid regime was accessed by keeping the 
lattice depth constant, for instance by ramping down the 
on-site energy using a Feshbach resonance. The Bragg 
signal as a function of the lattice depth is reported in 
Fig. [SI showing the appearance of the background signal 
as the superfluid regime is approached. 



IV. CONCLUSIONS 

We have discussed Bragg spectroscopy of ultracold 
bosonic atoms in an optical lattice, focussing on the sig- 
natures of the Mott-insulator and superfluid quantum 
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FIG. 5: (color online) Intensity of the scattered light (in arbitrary units) as a function of the Bragg angle G = Qxdo (in units 
of tt). The parameters are the same as in Fig. [2] and (a) Vo = S.IEr ( U/J « 17 ), (b) Vo = O.IEr ( U/J ^ 1). The black 
solid line corresponds to the numerical result, the blue dashed-dotted line to the analytical solution, the red dashed line to the 
numerical result obtained discarding the light-induced hopping term as in [24. i32i]. 




FIG. 6: (color online) Intensity of the scattered light (in ar- 
bitrary units) as a function of the Bragg angle O = Qxdo 
(in units of it) and of the lattice depth Vo (in units of hujR) 
(the corresponding values of the ratio U/J axe reported be- 
tween squared brackets). The other parameters are reported 
in Fig. [3 

state in the scattered photons. A full quantum theory 
for the atoms and photons dynamics and interactions has 
been developed, allowing us to identify the various con- 
tributions to the detected signals. We have characterized 
the Bragg scattering signal, for the parameters sweeping 
across the transition from the Mott-insulator to the su- 
perfluid quantum state. In particular, the contribution 
of light-induced hopping, arising from atomic recoil due 
to photon scattering, has been put into evidence. This 
term has been neglected in previous theoretical treat- 
ments [13, . In this work we have shown that its con- 
tribution can interfere with ordinary tunneling between 
sites thereby significantly affecting the spectroscopic sig- 
nal. Its effect is visible in the behavior of the height of the 



peaks in the spectrum as a function of the emission angle, 
and it has been singled out by comparing the spectrum 
evaluated when this effect is discarded. This effect can 
be revealed experimentally in large systems, according 
to the anal ytic al theory we develop by extending the one 
derived in |2^ . l4l| . and in small systems, as we observe 
by numerically evaluating the spectrum for a lattice of 7 
atoms. It is interesting to consider whether such proper- 
ties can be used as resources for photonic interfaces based 
on strongly-correlated atoms in optical lattices. 

This analysis has been made in the linear regime, as- 
suming a weak probe and far-off resonance both from 
the atoms and from the frequency of the lattice beam. 
Using instead Bragg beams at the same frequency as 
the optical lattice, wave-mixing effects are expected, 
as reported for instance in 0, [^, [5^. In addition, 
optical lattices have been discussed in the literature 
as a possible realization of p hotonic bandgap materi- 
als [il, m, H, [Hi, [13, m, ia [131 . An interesting ques- 
tion is how such photonic properties are modified when 
the many-body quantum state of the atoms is relevant to 
the atom-photon interactions dynamics. When the light 
is close to resonance with the atoms, hence in the dissipa- 
tive regime, the state of the atoms is significantly heated 
up. On the other hand, interesting photon-photon corre- 
lations could be observed, due to interference in multiple 
scattering by the atoms, see e.g. Ref. f6l| . 

We remark that, while monitoring the state of the gas 
by means of photons is attractive, on the other hand 
Bragg spectroscopy modifies the atomic system, as the 
recoil imparted by the scattered photon significantly per- 
turbs the state of the atomic gas. It would be desirable to 
identify schemes, such as quantum-non-demolition type 
of measurements [6^, [63j , which can allow one to measure 
the relevant quantities in a non-invasive way. This may 
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permit one to implement feedback mechanisms [36|, |6J] , 
which would allow one to prepare other nonclassical 
states of the atomic gas. 
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